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Abstract. We study charm production in ultra-relativistic heavy-ion collisions by using the 
Parton-Hadron-String Dynamics (PHSD) transport approach. The initial charm quarks are 
produced by the PYTHIA event generator tuned to fit the transverse momentum spectrum and 
rapidity distribution of charm quarks from Fixed-Order Next-to-Leading Logarithm (FONLL) 
calculations. The produced charm quarks scatter in the quark-gluon plasma (QGP) with the 
off-shell partons whose masses and widths are given by the Dynamical Quasi-Particle Model 
(DQPM), which reproduces the lattice QCD equation-of-state in thermal equilibrium. The 
relevant cross sections are calculated in a consistent way by employing the effective propagators 
and couplings from the DQPM. Close to the critical energy density of the phase transition, 
the charm quarks are hadronized into D mesons through coalescence and/or fragmentation. 
The hadronized D mesons then interact with the various hadrons in the hadronic phase with 
cross sections calculated in an effective lagrangian approach with heavy-quark spin symmetry. 
The nuclear modification factor Raa and the elliptic flow V2 of D° mesons from PHSD are 
compared with the experimental data from the STAR Collaboration for Au+Au collisions at 
y/SNN =200 GeV and to the ALICE data for Pb+Pb collisions at yj snn =2.76 TeV. We find that 
in the PHSD the energy loss of D mesons at high pr can be dominantly attributed to partonic 
scattering while the actual shape of Raa versus pr reflects the heavy-quark hadronization 
scenario, i.e. coalescence versus fragmentation. Also the hadronic rescattering is important for 
the Raa at low pr and enhances the D -meson elliptic flow V2- 


1. Introduction 

The Quantum Chromo Dynamics (QCD) [UEJG2I1] predicts that matter changes its phase at high 
temperature and density and bound (colorless) hadrons dissolve to interacting (colored) quarks 
and gluons, i.e. to the Quark-Gluon-Plasma (QGP). Such extreme conditions have existed in 
the early expansion of the universe and now can be realized in the laboratory by collisions of 
heavy-ions at ultra-relativistic energies. The study of the phase boundary and the properties of 
the QGP are the main goal of several present and future heavy-ion experiments at SPS (Super 
Proton Synchrotron), RHIC (Relativistic Heavy-Ion Collider), LHC (Large Hadron Collider) 
and the future FAIR (Facility for Antiproton and Ion Research) and NIC A (Nuclotron-based 
Ion Collider fAcility) [5|. Since the QGP is created only for a short time (of a couple of fm/c) it 
is quite challenging to study its properties and to find the most sensible probes. The advantage 
of ’hard probes’ such as mesons containing heavy quarks (charm and beauty) is, firstly, that 




due to the heavy masses they are dominantly produced in the very early stages of the reactions 
with large energy-momentum transfer, contrary to the light hadrons and electromagnetic probes. 
Secondly, they are not in an equilibrium with the surrounding matter due to smaller interaction 
cross sections relative to the light quarks and, thus, may provide an information on their creation 
mechanisms. Moreover, due to the hard scale, perturbative QCD (pQCD) should be applicable 
for the calculation of heavy quark production. As shown in Ref. [6], the FONLL calculations are 
in good agreement with the experimental observables on charm meson spectra in p+p collisions. 
This provides a solid reference frame for studying the heavy-meson production and their flow 
pattern in heavy-ion collisions. 

The first charm measurements at RHIC energies by the PHENIX[7] and STAR [8] 
collaborations were related to the single non-photonic electrons emitted from the decay of charm 
mesons. However, recently the STAR Collaboration measured directly the nuclear modification 
factor and the elliptic flow of D° mesons in Au+Au collisions at -^/snn =200 GeV |9J jTO] which 
allows for a straightforward comparison with the theoretical model calculations. It has been 
observed that the Raa and ^2 of charm mesons show a similar behavior as in case of light 
hadrons contrary to expectations from pQCD. Similar observations have been made at LHC 
energies, too mi¬ 
lt still remains a challenge for the theory to reproduce the experimental data and to explain 
simultaneously the large energy loss of charm quarks (Raa) and the strong collectivity (^ 2 ) (cf. 
e.g. [is ns nn ns]). Commonly, the interactions of charm quarks with the partonic medium 
are based on pQCD with massless light quarks and a fixed or running coupling. However, lattice 
QCD and quasiparticle approaches lead to the notion of massive degrees-of-freedom with finite 
spectral width m ■ In this case the pQCD cross sections are no longer meaningful. Moreover, 
the conclusions on the amount of suppression due to collisional energy loss by means of the 
elastic interactions of charm quarks with the QGP partons versus the radiative energy loss due 
to the emission of soft gluons (i.e. gluon bremsstrahlung) are still far from being robust. Also 
the influence of hadronization and especially hadronic rescattering is not yet settled, too. 

Our goal here is to study the charm dynamics based on a consistent microscopic transport 
approach for the charm production, hadronization and rescattering with the partonic and 
hadronic medium. In this study we will confront our calculations within the Parton-Hadron- 
String Dynamics (PHSD) approach to the experimental data on charm at RHIC and LHC 
energies and discuss the perspectives/problems of using the charm quarks for the tomography 
of the QGP. To achieve this goal we embed the heavy-quark physics in the existing PHSD 
transport approach m which incorporates explicit partonic degrees-of-freedom in terms of 
strongly interacting quasiparticles (quarks and gluons) in line with an equation-of-state from 
lattice QCD (1QCD) as well as dynamical hadronization and hadronic elastic and inelastic 
collisions in the final reaction phase. Since PHSD has been successfully applied to describe 
the final distribution of mesons (with light quark content) from lower SPS up to LHC energies 
BH ng EB !2Q], it provides a solid framework for the description of the creation, expansion 
and hadronization of the QGP as well as the hadronic expansion with which the heavy quarks 
interact either as quarks or as bound states such as D-mesons. We note in passing that in 
the hadronic phase PHSD merges with the familiar Hadron-String-Dynamics (HSD) approach 

mm- 

2. Description of charm production and dynamics in PHSD 

2.1. Charm production in p+p collisions 

Before studying the charm production in relativistic heavy-ion collisions, we discuss the charm 
production in p+p collisions at the top RHIC energy of -^/snn = 200 GeV. The charm production 
in p+p collisions also plays the role as a reference for the nuclear modification factor Raa in 
heavy-ion collisions. 
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Figure 1. (l.h.s.) Transverse momentum distributions (upper lines) of charm quarks in p+p 
collisions at y+NN =200 GeV from FONLL (dotted lines) and the tuned Pythia event generator 
(dot-dashed lines); in the lower lines the transverse momentum spectrum of D° mesons - which 
are fragmented from charm quarks with the contribution from D* decay included (solid line) 
- and that of D* + after scaling are compared with the experimental data from the STAR 
Collaboration [21]. (r.h.s.) The total cross section for charm production in p+p collisions 

(as parameterized in the PHSD) in comparison with the experimental data at various collision 
energies [211 125] , 


We use the Pythia event generator to produce charm and anticharm quarks in p+p collisions 
with the parameters PARP(91)=1.0 GeV/c and PARP(67)=1.0 as in Ref. [23]. The former 
parameter is the Gaussian width of the primordial transverse momentum of a parton which 
initiates a shower in hadrons and the latter the parton shower level parameter. We note that by 
an additional suppression of the transverse momenta of charm and anticharm quarks by 10 % 
and a suppression of rapidities by 16 % the transverse momentum spectrum as well as rapidity 
distribution of charm and anticharm quarks from the Pythia event generator are very similar 
to those from the FONLL calculations in p+p collisions at y+NN =200 GeV as shown in Fig. [T) 
(l.h.s.). Here the red dot-dashed lines are from the Pythia event generator (after tuning) and 
the blue dotted lines from the FONLL calculations, respectively. 

The produced charm and anticharm quarks hadronize by emitting soft gluons. The 
probabilities for a charm quark to hadronize into D°, D + , D*°, D * + , Df, and A c are, 
respectively, are specified in Ref. [23] together with the fragmentation function. The solid 
line in Fig. |T] (l.h.s.) shows the transverse momentum spectrum of D° mesons after charm quark 
fragmentation including the contribution from the decay of D*° and D* + . We can see that our 
results reproduce the experimental data from the STAR Collaboration [24] reasonably well. 

Since charm-quark production requires a high energy-momentum transfer, the number of 
produced charm quark pairs in relativistic heavy-ion collisions is proportional to the number 
of binary nucleon-nucleon collisions Nbi n - Since the probability to produce a charm quark pair 
depends on invariant energy and is less than that for primary hard collision in the PHSD, the 
binary nucleon-nucleon collisions producing charm quark pairs in PHSD are chosen by Monte 
Carlo from the ratio of the cross section for charm production in nucleon-nucleon collisions, 
(T^y/s), to the inelastic nucleon-nucleon cross section. The actual charm cross section cr^(y/s) 
(red solid line) is shown in Fig. 1 (r.h.s.) in comparison to the experimental data from [23, 25]. 







2.2. Heavy quark scatteruing in the QGP 

Quarks, antiquarks, and gluons are dressed in the QGP and have temperature-dependent 
effective masses and widths. In the DQPM, the mass and width of the light partons are given by 
thermal quantum-field theory assuming leading order diagrams but the strong coupling g 2 (T) 
is fitted to lattice data on energy and entropy densities, etc. [T5]. Note that a nonzero width of 
a parton reflects the off-shell nature as well as the strong interaction of the quasi-particle and 
finite life-time. 

The charm and anticharm quarks produced in early hard collisions interact with the dressed 
off-shell partons in the QGP. The cross sections for the charm quark scattering with massive off- 
shell partons are calculated considering the mass spectra of final state particles [2SU23- In the 
current study the charm quark mass is taken to be 1.5 GeV and its mass spectrum is neglected 
for simplicity. 

In the present study we consider only elastic scattering of charm quarks by light quarks and 
gluons. We do not consider yet the radiative processes which generate radiative energy loss 
because we expect that, due to the large gluon mass in the DQPM, the radiative processes are 
sub-dominant as compared to the collisional ones, especially for low charm-quark momenta (pt)- 
We expect the radiative energy loss to contribute at very high pt as accessible experimentally 
at the LHC [28] . 

We emphasize that the transport coefficient D s for charm quark diffusion - extracted from 
our microscopic calculations - and its agreement with the 1QCD results (within errors) and the 
corresponding D meson D s in hadronic medium validate our description for the coupling of 
charm with the QGP matter [29]. 

2.3. Hadronization of charm quarks 

Since the hot and dense matter created by a relativistic heavy-ion collision expands with time, 
the energy density of the matter decreases and the deconfined degrees-of-freedom hadronize to 
color neutral hadronic states. Once the local energy density in PHSD becomes lower than 0.5 
GeV/fm 3 , the partons are hadronized. First we look for all combinations of a charm quark and 
light antiquarks or of an anticharm quark and a light quark and calculate the probability for 
each combination to form a D or D* (or D S ,D*) meson. The probability for a quark and an 
anti-quark to form a meson is given by 
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where gM is the degeneracy of the meson M , and 
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with nii, r,; and k^ being the mass, position and momentum of the quark or antiquark i, 
respectively. The width parameter 5 is related to the root-mean-square radius of the meson 
produced through (r 2 ) = 3(m 2 + m 2 ) / (2(mi + m 2 ) 2 ) S 2 and thus determined by experiment (if 
available). However, we will use this radius as a free parameter and evaluate the coalescence for 
a D-meson radius of 0.5 and 0.9 fm (see below). For further details we refer the reader to [23]. 

Collecting all possible combinations of a charm or an anticharm quark with light quarks or 
antiquarks and calculating the coalescence probability for each combination from Eq. © , we 
obtain the probability for the charm or the anticharm quark to hadronize by coalescence in the 
actual space-time volume AtAxAyAz. Whether the charm or the anticharm quark is actually 
hadronized by coalescence is decided by Monte Carlo. Once the charm or the anticharm quark 
is decided to be hadronized by coalescence, then we find its partner again by Monte Carlo 






on the basis on the probability of each combination in the selected local ensemble. In case 
the charm or anticharm quark is decided not to hadronize by coalescence, it is hadronized by 
the fragmentation method as in p+p collisions (cf. Section 2.1). Since the hadronization by 
coalescence is absent in p+p collisions, it can be interpreted as a nuclear matter effect on the 
hadronization of charm and anticharm quarks. 

2-4■ D meson scattering in the hadron gas 

The D and D* mesons produced through coalescence or fragmentation interact with the 
surrounding hadrons in PHSD. The presence of several resonances close to threshold energies 
with dominant decay modes involving open-charm mesons and light hadrons suggests that the 
scattering cross sections of a D/D* off a meson or baryon, highly abundant in the post- 
hadronization medium, could manifest a non-trivial energy, isospin and flavor dependence. 
An example of these states is the broad scalar resonance D$ (2400), which decays into the 
pseudoscalar ground state D by emitting a pion in the s-wave (similarly to the heavy-quark 
spin partner Di(2420), decaying into D*n). Moreover, the similarity between the A(1405) and 
the A c (2595) has driven the attention to the fact that the latter could be playing the role of a 
sub-threshold resonance in the DN system, connected to the latter by coupled-channel dynamics. 

All these features have been addressed within several recent approaches based on hadronic 
effective models which incorporate chiral symmetry breaking in the light sector. The additional 
freedom stemming from the coupling to heavy-flavored mesons is constrained by imposing 
heavy-quark spin symmetry (HQSS). Whereas chiral symmetry fully determines the scattering 
amplitudes of Goldstone bosons with other hadrons at leading order in a model independent 
way, by means of HQSS the dynamics of the pseudoscalar and the vector mesons containing 
heavy quarks can be connected, since all kinds of spin interactions are suppressed in the limit 
of infinite quark masses [30] . For further details we refer the reader to Ref. [23j and the original 
literature cited therein. 


3. Comparison to experiment 

We now turn to actual nucleus-nucleus collisions at RHIC and LHC energies. Here the nuclear 
modification of D mesons is expressed in term of the ratio Raa which is defined as 


Raa {pt) = 


dNp U+Au 


/dp T 


^ dN^/dpr’ 


(3) 


where N^ u+Au and IV|) +P are, respectively, the number of D mesons produced in Au+Au 
collisions and that in p+p collisions, and is the number of binary nucleon-nucleon 

collisions in Au+Au (Pb+Pb) collisions for the centrality class considered. If the matter 
produced in relativistic heavy-ion collisions does not modify the D meson production and 
propagation, the numerator of Eq. ([3]) should be equal to the denominator. Therefore, an Raa 
smaller or larger than one implies that the nuclear matter suppresses or enhances D mesons, 
respectively. The elliptic anisotropy in the azimuthal angle is characterized by 


v 2 = (cos(2-0 - 2Tpp)} 


/ Pl-Py 

\P 2 x+p 2 y 



(4) 


where p x and p y are the x and y components of the particle momenta and \kpp is the azimuth 
of the reaction plane while the brackets denote averaging over particles and events. The v 2 
coefficient can be considered as a function of centrality, rapidity y and/or transverse momentum 
Pt- We note that the reaction plane in PHSD is given by the (x — z) plane with the z-axis in 
the beam direction. 






Figure 2. (l.h.s.) The Raa of D° 

mesons including partonic scattering with 
(dashed and solid lines) and without hadronic 
scattering (dotted and dot-dashed lines) in 
Au+Au collisions at y+NN =200 GeV for 
a D meson radius of 0.5 fm and of 0.9 fm. 
The experimental data are from the STAR 
Collaboration II. (r.h.s.) The elliptic 
flow V 2 of D° mesons including partonic 
scattering with (dashed and solid lines) 
and without hadronic scattering (dotted and 
dot-dashed lines) in Au+Au collisions at 
y+nn =200 GeV for a D meson radius of 
0.5 fm and of 0.9 fm. 


Fig. [2] shows the Raa (l.h.s.) and the elliptic flow V 2 of D° mesons (r.h.s.) with and without 
hadronic scattering. We can see that the hadronic scattering plays an important role both in 
Raa and the elliptc flow V 2 - It shifts the peak of Raa to higher transverse momentum especially 
in central collisions and enhances the elliptic flow of final D mesons. 

In Fig. [3] we present our preliminary results for the LHC energies. The left part of the 
figure shows the PHSD results for the Raa of D° + D + + D* + mesons at midrapidity for 20% 
central Pb+Pb collisions at y+NN =2.76 TeV calculated without (red line) and with (blue line) 
shadowing m in comparison to the experimental data from the ALICE Collaboration [32] • As 
seen the shadowing effect reduces the peak at low px- which is still higher than the experimental 
data. The suppression of the high px tail is comparable with the corresponding results for RHIC 
energies and approximately in line with the data. We stress again that we did not account for 
radiative energy loss in the present calculations which might play a role at higher px [28]. 

The right part of figure [3] shows the corresponding V 2 for 20% central Pb+Pb collisions for 
the scenario including shadowing. We note that the t> 2 -fl° w for central (0-20%) collisions is 
lower then that for semi-central (30-50%) measured by the ALICE Collaboration [333, however, 
still quite substantial (up to 10%) and comparable with the V 2 of light hadrons for the same 
centrality. 

4. Summary 

We have studied charm production in relativistic heavy-ion collisions by using the Parton- 
Hadron-String Dynamics (PHSD) approach [18] where the initial charm quark pairs are produced 


































Figure 3. (l.h.s.) The PHSD results for the Raa of D° + D + + D* + mesons at midrapidity 

for 0-20% central Pb+Pb collisions at y / SNN =2.6 TeV calculated without (dashed blue line) 
and with (solid red line) shadowing in comparison to the experimental data from the ALICE 
Collaboration |32j. (r.h.s.) The PHSD results for the V 2 versus px of D° + D + + D* + mesons 
for 20% central Pb+Pb collisions at ^/snn =2.6 TeV including the shadowing effect. 


in binary nucleon-nucleon collisions from the PYTHIA event generator taking into account the 
smearing of the collision energy due to the Fermi motion of nucleons in the initial nuclei. The 
produced charm and anticharm quarks interact with the dressed quarks and gluons in the QGP 
which are described by the Dynamical Quasi-Particle Model m in PHSD. The interactions of 
the charm quarks with the QGP partons have been evaluated with the DQPM propagators and 
couplings consistently [26]. We mention that when extracting the spatial diffusion coefficient 
D s from our cross sections as a function of the temperature we observe a minimum of D s close 
to T c which is in line with lattice data above T c and hadronic many-body calculations below T c 
(cf. Ref. [29]). 

We have found that the interaction of the charm quarks with the dynamical partons of 
the QGP softens the px spectrum of charm and anticharm quarks but does not lead to a full 
thermalization for transverse momenta px > 2 GeV/c. The charm and anticharm quarks, 
furthermore, are hadronized to D mesons either through the coalescence with a light quark 
or antiquark or through the fragmentation by emitting soft ’perturbative’ gluons. Since the 
hadronization through coalescence is absent in p+p collisions, it can be interpreted as a nuclear 
matter effect on the D meson production in relativistic heavy-ion collisions. In the coalescence 
mechanism the charm or anticharm quark gains momentum by fusing with a light quark 
or antiquark while it looses momentum in the fragmentation process (as in p+p reactions). 
This partly contributes to the large Raa of D mesons between 1 and 2 GeV/c of transverse 
momentum. Finally, the formed D mesons interact with hadrons by using the cross sections 
calculated in an effective lagrangian approach with heavy quark spin-symmetry [30], which is 
state-of-the art. We have found that the contribution from hadronic scattering both to the Raa 
and to the elliptic flow of D mesons is appreciable, especially in central collisions, and produces 
additional elliptic flow V 2 [23]. 

Since the PHSD results reproduce the experimental data from the STAR and ALICE 
Collaborations without radiative energy loss in the px range considered, we conclude that 
collisional energy loss is dominant at least up to px = 6 (15) GeV/c in relativistic heavy- 
ion collisions at RHIC (LHC). In our approach this is essentially due to the infrared enhanced 
coupling a s (T ) in the DQPM leading to large scattering cross sections of charm quarks with 
partons at temperatures close to T c . 
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